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Abstract 

The gravitational wave (GW) memory from a radiating and decelerating point mass is studied in detail. 
It is found that for isotropic photon emission the memory generated from the photons is essentially the 
same with the memory from the point mass that radiated the photons so that it is anti-beamed. On the 
other hand, for anisotropic emission the memory from the photons may have a non-vanishing amplitude 
even if it is seen with small viewing angles. In the decelerating phases of gamma-ray burst (GRB) jets the 
kinetic energy of the jet is converted into the energy of gamma-ray photons. Then it would be possible 
to observe a variation in the GW memory associated with GRB jets on the timescale of the gamma-ray 
emission if the emission is partially anisotropic. Such an anisotropy in the gamma-ray emission has been 
suggested by the polarizations detected in recent observations of GRBs. The GW memory from GRB jets 
would provide clues to clarifying the geometry of the jets and the emission mechanism in GRBs. Thus it 
will be an interesting target for the next generation detectors of the GWs. 
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1. Introduction 

Gamma-ray bursts (GRBs) are considered to be rela- 
tivistic jets with Lorentz factors 7 ~ 10^ — 10"^. The en- 
ergy radiated in gamma-rays for a single burst amounts to 
~ 10^° - 10^^ erg (Frail et al. 2001), which may be sup- 
plied by kinetic energy of the relativistic jets (e.g.,Piran 
1999). The origin of such energetic jets currently de- 
bated includes the coalescence of double neutron stars 
(NSs) in a binary as well as the collapsar, which is the 
collapse of a rotating massive star ending up with a sys- 
tem composed of a black hole (BH) and an accreting torus 
(Woosley 1993). It has long been considered that NS bi- 
naries and collapsing massive stars are the primary tar- 
gets for the gravitational wave (GW) detection for mis- 
sions such as LIGO, TAMA, Virgo, and LISA. In fact, 
some experiments have already begun to put constraints 
on the rates of such astrophysical phenomena(Akutsu et 
al. 2006; Harry et al. 2010). Soon after the discovery of 
the supernova-GRB association (Galama et al. 1998), a 
hydrodynamical simulation for the collapsar model was 
carried out by MacFadyen & Woosley (1999) and later 
it has been studied intensively by a number of numerical 
simulations (Nagataki 2009; Sekiguchi & Shibata 2011; 
Ott et al. 2011 to list only a few of the recent works). It 
is now established that the collapsar scenario is the most 
plausible model that wc now know for long GRBs. 

Then GRBs turn out to be major astrophysical sources 
of GWs. The mechanism of GW emission from GRBs 
is quite similar to that of the GWs from supcrnovae, 
which is mainly caused by asphericity in the matter mo- 



tion in a bouncing core of collapsing massive stars and 
the anisotropic neutrino emission that emerges from the 
core (Miiller & Janka 1997; Dimmelmeier et al. 2007). In 
the collapsar scenario of GRBs, the neutrinos are emitted 
from an accretion torus and the total energy of neutri- 
nos are estimated to be - 10^^ _ iQ^t q^^^^ q^^ 
signals have been studied analytically (Hiramatsu et al. 
2005; Suwa & Murase 2009) and by special relativistic 
magnetohydrodynamics simulations (Kotake et al. 2012). 
The GW from anisotropic neutrino emissions comes up 
as a burst with memory, that is, a sharp rise followed by 
a steady value of the metric perturbation (Braginsky & 
Thorne 1987). In fact the GW memory has been con- 
sidered for a variety of sources such as relativistic fluids 
(Segalis & Ori 2001; Sago et al. 2004) and neutrinos ( 
Hiramatsu et al. 2005; Suwa & Murase 2009) in the con- 
text of the GW from GRBs. The maximum amplitude of 
the GW memory produced by a relativistic point source 
with energy E is given by 



A/i= — — 



1.1 X 10" 



E 



lOSiergy VlMpc 



(1) 



where d is the distance to the source. Here we mean 
by a 'point source' a point mass or a single ray of pho- 
tons/neutrinos, where the point mass is an approximation 
to treat a fluid element with a non-zero mass density. 

Segalis & Ori (2001) studied the GW memory from a 
relativistic jet of GRBs. They derived the angular de- 
pendence of the GW amplitude for the case of a point 
mass and found that the amplitude becomes quite small 
in the direction of jet propagation (anti-beaming effect). 
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Another important finding is that the GW memory is po- 
larized in the direction from the source to the line of sight 
on the transverse plane ( in the TT gauge ). Sago et al. 
(2004) analyzed the GW memory from the accelerating 
phases of the GRB jets. They find that the finite size of 
the opening half-angle /S.0 of the jet leads to the reduction 
of the GW memory for viewing angles 9^ smaller than A6', 
which means that the GW memory is hard to detect simul- 
taneously with GRBs. The total energy of neutrinos is two 
orders of magnitude larger than in the case of relativistic 
jets of GRBs. However, the anisotropy in the neutrino 
emission is weak compared to GRB jets, resulting in the 
same consequence that it is not likely that we could ob- 
serve the GW memory at the same time with GRB events 
(Hiramatsu et al. 2005; Suwa & Murase 2009). 

The amplitude of the GW memory expected for the rel- 
ativistic jet of GRBs is much smaller than that expected 
for neutrinos from GRBs. However, a possible detection of 
the GWs characteristic for the jets would support the gen- 
eral view that GRBs are associated with relativistic jets 
and might resolve the issues on the physical mechanism 
of the gamma-ray emission. Sago et al. (2004) investi- 
gated the GW memory from relativistic jets based on the 
'unified model' of GRBs , in which a bunch of multiple 
sub-jets launched from a central engine are assumed to 
be seen as GRBs, X-ray flushes (XRFs), and X-ray rich 
GRBs, depending on the viewing angle (Yamazaki et al. 
2004). They also pointed out that GRB jets may be GW 
sources in deci herz frequency bands, reflecting the char- 
acteristic time scale of the central engine's activity, and 
thus will be suitable targets for DECIGO and BBO (Seto, 
Kawamura, & Nakamura 2001; Sago et al. 2004). 

It is thought that GRB phenomena detected by GW 
memory will be seen from ofF-axis with Oy > because 
of the anti-beaming effect. Thus, such phenomena should 
actually be observed as XRFs rather than GRBs, or com- 
pletely be missed, in electromagnetic observations. In the 
GRB phenomena, the relativistic jet will lose its kinetic 
energy and will be decelerated because of the gamma-ray 
emission. Then part of the memory carried initially by 
the jet will be transformed into the memory produced 
by electromagnetic radiation. Although the anti-beaming 
effect does not exist for the memory generated by a sin- 
gle ray of photons, photons are emitted toward a small 
but finite solid angle with a typical opening half-angle of 



celerating point mass. In §3, we apply the result of §2 
to realistic cases corresponding to GRB phenomena. The 
finiteness of the opening half-angle of the jet is now taken 
into consideration. Then our analyses will be applied to 
the internal shock model of GRBs. We also present cal- 
culated waveforms for a specific model that has appeared 
in the previous literature. Finally, the summary and dis- 
cussion including the detectability are given in §4. 

2. Gravitational Wave Memory from a Radiating 
Point Mass 

We assume a flat background space-time with the 
Minkowski metric 77^,^ = (— l,+l,-fl,-|-l) and restrict 
ourselves to dealing with the small perturbation h^^ = 
gfiv ~ Vfj.i' of the metric g^jy. We set G = 1 hereafter. 
Then the linearlized Einstein equation reads 



(2) 



under the Lorentz gauge condition hf^"^^ = 0, where ft.^^ is 
defined as 

7;V,^^h^, (3) 



and Tfj_i, is the energy-momentum tensor. For a point mass 
m with a world-line x^{t) parametrized by the proper 
time T, T'^'^ is given by 



T^"'{x)= / mu^'{T)u''{T)5^'^\x~x{T))dT. 



(4) 



where u'^{t) = ~ 7(1, /3n) and /?, n are the velocity, 
the unit vector proportional to the velocity of the point 
mass, respectively. 

Following the notation in Segalis & Ori (2001), the re- 
tarded solution of equation (2) is written as 



U^'{Tr)u''{Tr) 



(5) 



-Ux{Tr) ■ {x-x{Tr))^' 

where Tr is the retarded time defined by the conditions 

{x - a;(T^))^ • (x - x{t,.))^ =0, x° - ^^(t^) > 0(6) 

The transverse-traceless (TT) part of hij defined by 
hf^ = PikhkiPij - \PijhkiPku where = 5ij - n[n'^ and 
n'^ is i—ih component of the unit vector from the source 
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so that the phase cancellation of the amplitude still to the observer, has been derived to be 



occurs. However, if the emission is partially anisotropic, 
we might be able to see a variation in the GW memory 
from GRB jets for a moderate range of viewing angles 

^ 7^^ ^ on the time scale of gamma-ray emission. 
In fact, the polarizations detected in recent observations 
of GRBs suggest that the gamma-ray emission is partially 
anisotropic (Steele et al. 2009; Yonetoku et al. 2011). 

In this paper, we will present the results of a detailed 
study on the change of the GW memory in the decelerat- 
ing relativistic shocks of GRBs. The paper is organized 
as follows. In §2 we briefly describe the formulation for 
calculating the gravitational waveform and focus on the 
change of GW memory generated by a radiating and de- 



lTT 



"-yy 
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27771/^2 



d 1 — /ScosS 
27TO/32 sin^ 6 



cos2(/). 



d 1 — /3cos6' 



sm2(f>, 



(7) 



in Segalis & Ori (2001). Here we choose z axis to be the 
line of sight, d is the distance to the source and 6,4> are 
the polar and azimuthal angle in spherical coordinates, re- 
spectively, that indicate the direction of the jet. It should 
be noted that /3 and 7 are evaluated at the retarded time. 
It is found that the 9 dependent part of and , which 
we denote h{9), has asymptotic forms 



No. 



HO) 



sm 
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(8) 



l-ZScosS 7-2 + 612 

for 7 >> 1 and 9 << 1 and 

h{9)c^ l + cosO (9) 

for 9 large compared to 7~^. This angular dependence 
shows that the GW memory is diminished in the forward 
direction within < 7"-'^ , which is called as " anti-beaming" 
(Segalis & Ori 2001; Sago et al. 2004). 

Now we turn to the case of photons. Since there is no 
rest frame for a zero-mass particle, the world-line is a null 
geodesic, being parametrized as x'^{x°) by the coordinate 
time a;°. Then the energy momentum tensor is expressed 
similarly as 



T^^-'ix) = En^'{x°)n''{x°)S^^^ (x - x(a;")) 
En''{t)n''{t)5'-^\x-x{t))dt 



(10) 



where a;°(i) = t, E is the energy of the photon, n'^ = (l,n) 
is the unit four- vector that is proportional to the photon's 
four momentum. 

Comparing equations (4) and (10), we notice that the 
retarded solution for a zero-mass particle is obtained 
straightforwardly from the one for a point mass with a 
simple replacement 7m ^ E, /3 — ?> 1, ^ n^^ . Then we 
have the TT part of the metric perturbation as follows. 



2E 



(1- 



cos 6*) cos 20, 



) sin 20. 



(11) 



This expression is consistent with the formula given in 
Miiller & Janka (1997), which was derived for the GW 
memory from supernova neutrinos based on the work of 
Epstein (1978). It has already been suggested that the 
above replacement seems to be valid in Sago et al. (2004). 

The GW memory from photons ( equation 11) has the 
same angular dependence with the memory from a point 
mass for 9 > (equation 9). However, for small 9, the 
anti-beaming is not seen for photons unlike the case of a 
point mass. Then, we might think of the idea that we 
would observe an increase in the GW memory, for small 
viewing angles, under the conversion of energy from a 
point mass into photons. However, photons are emitted 
from the point mass into a finite solid angle with a typi- 
cal opening half-angle of ~ 7~^, which results in the same 
effect as if the jet has a finite opening half-angle. That 
is, the amplitude is likely to be canceled for small 9 if the 
emission would be axisymmetric. 

To calculate the GW memory, it is useful to define an- 
other frame {x' ,y' ,z'), where j/'-axis coincides with y-axis 
and the a;'(z')-axis is obtained by rotating the a;(z)-axis 
around the y axis by an angle 9y so that z' axis corre- 
sponds to the forward direction of the point mass or the 
symmetry axis of the distribution of the photon emission 
(figure 1). Thus 9y is equal to the viewing angle of the jet. 
In changing coordinates of both frames into the spherical 
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Fig. 1. Two coordinate frames (x,y,z) and {x' ,y' ,z'). The 
point mass (jet) moves in tlie z' direction and tlie 2 axis cor- 
responds to the line of sight. The viewing angle is indicated 

by e^j. 

coordinates, the following relations hold among 9,(f) and 



sm f cos (p = — sm 
sin 9 sin cb = sin 9' sin 



cost 
/ 



cosf,, smf cose 



cost 



cos t/„ cos 



9' + sin 9,, sin 9' cos c 



(12) 



Then, the GW memory should be averaged over the 
angular distribution of photons as follows. 



h[9y) = h+{9y) + ihy, (9,,) = 4 < 2(1 + cos6i)e^''^ > 

a 



2E 



sin9' d9' d(j)' f{9' ,4,'){l + cos9)e^^'l' (13) 



where 9 and are related to 9' and 0' by equation (12), 
E is the total energy of the photons, and f{9',(l>') is the 
distribution function. The fraction of the energy of pho- 
tons emitted into a solid angle d^l' = sm9'd9'd(j>' is given 
by f{9',(f>')dfl'. Figure 2 shows a schematic view of the 
cone of emitted photons. 

We evaluate the GW memory from photons for sev- 
eral cases of different angular distributions of the pho- 
ton emission. First we assume an isotropic emission in 
the rest frame of the point mass for simplicity. Some 
GRBs have shown relatively high degrees of polarization 
in their early optical emission (e.g., Steele et al. 2009) and 
in gamma-ray emission itself (e.g.,Yonetoku et al. 2011). 
Such high degrees of polarizations may suggest the exis- 
tence of ordered magnetic fields in the relativistic shock. 
If the gamma-rays are caused by synchrotron radiation 
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Fig. 2. Schematic view of the cone of emitted photons 

in ordered transverse magnetic fields (Granot & Konigl 
2003), the emission would tend to be focused into solid 
angles with a limited range of azimuthal angles so that 
the angular distribution will deviate from the axisymme- 
try. To model such cases in a simplest way, we study 
two cases. One is a maximally anisotropic case where the 
photon emission is focused in the direction of a particu- 
lar azimuthal angle. The other is the case of synchrotron 
emission in the presence of ordered magnetic fields which 
has been considered in the modeling of the polarization of 
the gamma-ray emission(Granot & Konigl 2003; Toma et 
al. 2009). 

2.1. Case oj Isotropic Emission 

If the photon emission is isotropic in the rest frame of 
the point mass, the angular distribution function is given 
by 

f{e',(j>')dn' — \ — —dn' (14) 

^ ' 47r74(l-/?cos6'')^ 

and /?, 7 is the velocity and the Lorentz factor of the 
point mass, respectively (Rybicki k. Lightman 1979). For 
7 > > 1 the emission is beamed within small angles of or- 
der 6' ^ 7^^ so that the amplitude depends only on 70„. 
Figure 3 shows two polarization components /i+ and 
of the GW memory from photons for isotropic emission 
generated by a point mass. The normalized amplitude 
h/{E/d) is shown as a function of ^6^ (solid line). The 
/ix component vanishes owing to the axisymmctry. We 
also plotted the memory expected from a point mass with 
the same energy 7m = E (dotted line). Since /i+ takes 
almost the same value with the memory from the point 
mass, it is hard to distinguish them in Figure 3. 

2.2. Case of Maximally Anisotropic Emission 

Next we consider the second case in which the angular 
distribution is given by 

1 1 



f{e',cty')dn' = 



-M'-Mdn\ (15) 



2 74(l-/3cos6l')^ 

where P and 7 is the velocity and the Lorentz factor of 
the point mass that radiates the photons. 

For 0„ = the GW memory from photons with the an- 
gular distribution of equation (15) can be evaluated ana- 
lytically as 



2E{1 + I3) 



(16) 
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Fig. 3. GW memory from photons with isotropic angular 
distribution of equation (14) radiated by a point mass with 
7 = 100. The amplitudes normalized hy E/d are shown as a 
function of the viewing angle 6„ (solid lines). The hx com- 
ponent vanishes owing to the symmetry. The memory from 
a point mass with the same energy 7m = E is also plotted 
(dotted line), which takes almost the same values with h-^- so 
that it is hard to distinguish them. 

For large 7 ~ 10^ — 10'^ the amplitude of the GW memory 
takes a maximum value and the polarization is simply 
determined by the angle 0o. 

For 6y 7^ we calculated the amplitudes averaged with 
equation (15) numerically. Figure 4 shows the two com- 
ponents of polarization and hx calculated for (f>o = 
O7 f I f I f a function of 76*^. Figure 5 shows the same 
amplitudes but for 0o = f j f ""i f'^'i I'"'- The amplitudes 
are normalized by E/d again. The GW memory for a 
point mass with jm = E is also shown with a dotted line 
for comparison. As 9y increases, becomes a dominant 
component, approaching to the memory for a point mass 
asymptotically. 

2.3. Case of Synchrotron Emission with Ordered 
Magnetic Fields 

Here we consider synchrotron emission from electrons 
as a more physically-motivated model. We assume that 
an ordered magnetic field B = (i? cos 0b, i? sin 0^,0) is 
present ( Figure 6 ), the emitted photons have a simple 
power-law energy spectrum Pi, oc and that the elec- 
trons have an isotropic pitch angle distribution (Granot 
& Konigl 2003). The photons emitted in the direction 
n= (sin6'' cos 0', sin0' sint/)', cos 0') have a pitch angle ^ 
B • n 

such that cos^ = — '^tiere represents the direc- 
tion of the photons in the comoving frame, 



lie = (sin cos (/)', sin 0c sin (/)', cos 6'c), 



(17) 
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Fig. 4. GW memory from photons with anisotropic angular distribution of equation (15) normalized by E/d (case of maximally 
anisotropic emission). Two polarization components /i+ and hx (solid lines) are shown for </>o = (top left), ^ (top right), -j 
(bottom left), and ^tt (bottom right) as a function of 76'^. The amplitude from a point mass with the same energy 7m = E is shown 
for comparison (dotted line). 




Fig. 5. The same with Fig. 4 but for 4>o = ^ (top left), ^tt (top right), ^tt (bottom left), and ^tt (bottom right). GW memory from 
photons with anisotropic angular distribution of equation (15) normalized by E/d (case of maximally anisotropic emission). The 
amplitude from a point mass with the same energy ■ym = E is shown for comparison (dotted line) . 
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Fig. 6. Geometry considered for the synchrotron emission 
from ordered magnetic fields. 



and 9c is given by cosOc = 



1 - {^B' f 



cosO' -P 
1-/3 cos 61'' 



which reduces to 



for 7 >> 1 and 9' «l. We note that 



1 + {i9'f 

the magnetic field in the comoving frame is given by 7B as 
a result of Lorentz transformation. A simple calculation 
yields 



COS^: 



: smfc cose, 
■ sv[\9r cosf 



cos (pB - 



sm^cSm^ su\<pB 



(18) 



Since the synchrotron power is proportional to 
(sin^)°+^ (Rybicki & Lightman 1979), the angular dis- 
tribution is given by 



f{9',4>')dVl' ^K{sin£,) 



a+l 



1 



(19) 



7^(1 cos 6*')^ 

where isT is a normalization constant. We choose a — \ 

from a range of values considered in the literature (Granot 

& Konigl 2003; Lazzati 2006; Toma et al. 2009). For a = 1 

3 

the distribution function has a simple form and K = — . 

Stt 

The GW memory depends on (j)B as well as 9y so that we 
denote the amplitude as h{9v,(j)B) in this case. For 9^ = 
the integration in equation (13) can be done analytically, 
where we find 



MO>b) 



1 2E{1 + P) .^^^ 
4' d 



(20) 



We see that 4>o dependence of h{0) in equation (16) is 



reproduced by replacing (j)B with < 



in equation (20). 



For 9v ^ we calculated the amplitudes averaged with 
equation (19) and a = 1 numerically. 

Figure 7 shows the normalized amplitudes h{9v,(t>B) = 
h+{9y,(j)B) + ihxi9v,4'B) calculated for 0b =0, |,|,|7r 



as a function of ^9y. Figure 8 shows the same for (pB ~ 
^, |7r, |7r, Itt. The GW memory for a point mass with 
7m = ii^ is also shown with a dotted line. The overall be- 
havior is quite similar to the case of maximally anisotropic 
emission except that the amplitude for 0^, = is reduced 
to one fourth of the maximum amplitude. 

Now we consider the variation in the GW memory as 
a result of photon emission from a point mass. If a point 
mass m with a Lorentz factor 7 ( velocity /?) radiates 
photons with power P during a short time interval At, 
the change in the total GW memory Ah = Ah^ + iAh^ 
is given by 



Ah: 



2j'm P''^sm'^9y PAt 
d 1 — P' cos 9i, d 
27m I3'^s\u^9y 



<2(1 



- cosf 



)e2# > 



d 1 — /3cos6 



(21) 



where 7'(/3') is the Lorentz factor (velocity) of the point 
mass after time At. Since the conservation of energy im- 
plies 7'm = 7m — PAt, wc find that the following relation 
holds to first order in At, 



27m , / fl^ sin^ 9^ 
Ah=——A' 



d \1 — /3cos6 
2PAt (2-cos6'.„)sin2 9y 



2l3'^sm'^9y 
cos 9„ 



d -f^{l~pcos9y)^ 

P^t ( , , 2sin26'„ , , , 

< 2(1 + cos0)e2^^ > ^ ,(22) 



1 



Jcos0« 



where we assume 7,7' >> 1 and retain the terms in lead- 
ing orders of 7~^. In equation (22) the first term is the 
variation in the memory of the point mass due to broad- 
ening of the hole of anti-beaming, while the second term 
is the change in the memory ascribed to the emission of 
photons. The size of the terms in the parenthesis of equa- 
tion (22) can be easily read from figures 3, 4, 5, 7, and 8 
as the difference between /i+ (solid line) and the memory 
from a point mass (dotted line) as well as the value of 
ft-x- For the case of isotropic emission, the second term 
is negligibly small compared to the first term ( Figure 3). 
Thus the memory decreases monotonically as the point 
mass slows down, even if the contribution from photons is 
added. However, for the case of anisotropic emission above 
considered (equations 15 and 19), the second term in equa- 
tion (22) makes a dominant contribution to the change in 
the memory compared to the first term for small 7^^ . For 
large ^9^ the first term in equation (22) is smaller than 
the second term in the parenthesis by a factor of 7^ and 
thus it is negligible again. As a consequence, for ^9^ < 1 
and even for 9y = the GW memory varies by a significant 
fraction of its maximum value on the conversion of energy 
from a point mass into photons. 
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Fig. 7. GW memory from photons with anisotropic angular distribution of equation (19) normalized by Ejd (case of synchrotron 
emission with ordered magnetic fields). Two polarization components and hx (solid lines) are shown for ijio = (top left), ^ 
(top right), J (bottom left), and |7r (bottom right) as a function of 701, . The memory from a point mass with the same energy 
7m = E is shown for comparison (dotted line) . 
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Fig. 8. The same with Fig. 7 but for 4>b = ^ (top left), |7r (top right), ^tt (bottom left), and ^tt (bottom right). GW memory 
from photons with anisotropic angular distribution of equation (19) normalized by E/d (case of synchrotron emission with ordered 
magnetic fields) . The amplitude from a point mass with the same energy 7m = E la shown for comparison (dotted line) . 
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3. Gravitational Wave Memory from the 
Decelerating Phase of GRB Jets 

Here we estimate the overall behavior of time variation 
in the GW memory from GRB jets by applying our re- 
sults in §2 to specific models of GRB phenomena. It has 
been shown that the temporal structure of GRBs are nat- 
urally explained by the internal shock model, in which 
gamma-rays are radiated in an inhomogcncous relativis- 
tic wind possibly generated by a variable central engine 
(Piran 1999). Such a wind has been modeled by multiple 
shells moving in procession with various Lorentz factors. 
It is assumed that if a rapid shell catches up a slower 
shell ahead the two shells would collide and merge con- 
verting a fraction of kinetic energy into thermal energy. 
The thermal energy released will be radiated on the spot. 
We follow the formulation by Kobayashi et al. (1997) to 
model the internal shock. Here we use the internal shock 
model only to reproduce the time variability of typical 
GRB light curves. 

We assume that a rapid shell with mass and Lorentz 
factor 7r collides with a slow shell with mass and 
Lorentz factor 7^ to form a merged shell with mass 
nir + ms and Lorentz factor 7^. The conservation of en- 
ergy and momentum leads to 

rur^r + msjs = {rrir + ms + 6)7^, 



rur \Jil-\ + ms - 1 = ("V + ms + e) V7|~T(23) 

where e is the internal energy being released in the rest 
frame of the merged shell. For large 7r,7s >> 1, 7m is 
approximately given by 



It, 



1 rrir^r - 




TTlrhr ~ 





(24) 



and then the energy converted into radiation E is esti- 
mated to be 



E = 7„ie = ™r(7r - 7m) + ^-3(75 - 7m)- (25) 

At the time of collision the forward and reverse shocks 
arise. The Lorentz factors of the forward and reverse 
shocks, 7/s and 7^3, arc given by 



If 8 ^ In 



Jrs — Jm 



respectively (Sari & Piran 1995). The time in which 
gamma-rays are emitted is approximately given by the 
time in which the reverse shock traverses the rapid shell, 

where Ir is the width of the rapid shell and the /3rs is 
the velocity of the reverse shock. We note that the emit- 
ting region moves at a speed of /3m and that the observed 




timescale is reduced by a factor of 1 — /3m cos 6y, which is 
^1/27^ for 0„=O. 

The width of the merged shell Im is simply given by 
widths of the rapid and slow shells lr-,ls as 



^m, ^. 



13 fs ~ , Pm— Pr 



(27) 



Pfs ~ Ps Pr~ Prs 

In calculating the GW memory, we assume that for each 
collision the energy of photons E (equation 25) is radiated 
in time At (equation 26) so that the emission power is 
given by P = E / /S.t. Applying equation (21) in the last 
section, the change in the GW memory for a single colli- 
sion is evaluated as 

A/l(6l„ , 0s ) = A/lphotons [Ov Ab) + A/ljet {9^), (28) 

where 



A/lphotons(^'t),0B) 

_ mr{"fr - 7m) + rns{"fs - -1m) 

d 

for photons' contribution and 



<2(l + cos6l)e2"^ >(29) 



A/i 



2TOoSin^ 9, 



2mr sin^f 



ImPn 



IrPl 



d \ 1 — 13m COS 9„ 

ImPm IsPs 



1 ~ (3r COS i 



(30) 



d V 1 ~ Pm COS 9y 1— Ps COS 9y 

for the fluid contribution in the shells. We assume that 
this change in the memory occurs during time At at a 
constant rate. We study the case of synchrotron emission 
with ordered magnetic fields for A/iphotons- We always 
have a real and negative A/ijot, while A/iphotons may be 
complex reflecting the non-axisymmctric distribution of 
the photon emission. 

Now we describe the method for simulating the collision 
of the multiple shells in the internal shock. We take N 
shells with the same mass m and width I (i.e., the same 
density) being placed with an equal separation L at time 
t~0. To each shell a Lorentz factor, uniformly distributed 
between 7,„i,i and 7max, is assigned at the outset. Given 
the Lorentz factor 7^, the velocity and the position 
Xi = —{L + l)i for i-th shell ( i = 1, • • • ,iV ) as the initial 
conditions, we follow the time evolution of the position 
of all the shells until the first collision occurs. At the 
collision two shells are merged and again the movement 
of the shells is followed till the next collision. The same 
procedure is repeated until either all shells will be merged 
into a single large shell or we get the configuration where 
the shells are marching with decreasing velocities from the 
head to the end. We used parameters N = 100, 7min = 
10^,7max = 10'^, and i = 1,Z = 0.1. The overafl evolution 
of the GW memory witli time and the gamma-ray light 



curve arc determined by the ratio L/l, iV, 7n 



7n 



and 



Since the GW memory is polarized in the direction from 
the source to the line of sight on the transverse plane, we 



No. ] Gravitational Wave Memory from the Relativistic Jet of GRBs 




^ y' 



X 



Fig. 9. A uniform jet with an opening half angle A9 is seen 
from the top of the jet. Only a small portion of the jet around 
the line of sight of an angular size (shaded circle) con- 
tributes to the change in the GW memory. 
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been considered to explain a large degree of polarization 
(Granot & Konigl 2003; Lazzati 2006; Toma et al. 2009). 
If the angular scale of coherent magnetic field is small 
compared to the beaming angle 7"^ or the magnetic fields 
are randomly oriented, there is no preferred azimuthal an- 
gles into which the gamma-ray emission is focused so that 
the emission ought to be essentially isotropic. Then the 
change in the GW memory on each collision should be 
quite small so that we would not observe a variation in 
the memory from the decelerating phases of GRB jets. 

Instead if the shell contains a number of patches of 
locally ordered magnetic fields we would have a non- 
vanishing variation in the memory for some viewing an- 
gles. The case in which we would expect the largest change 
in GW memory is when an ordered transverse magnetic 
field exists in the relativistic shell. We consider a uniform 
jet with an opening half angle A9 illustrated in Figure 9. 
The line of sight ( z axis ) is located at a point tilted from 
z' axis by an angle 9y toward x' direction. We assume 
that the magnetic field has only a transverse component 
and that the field line is parallel to the meridional plane 
at (j)' ~ a. Then the unit vector B = B/|B| at a point 
[9' ,<f)') is expressed as 



need to account for the finiteness of the conical opening 
half-angle of the GRB jets to estimate the net memory. 
The typical size of the opening half-angle has been esti- 
mated to be ~ 0.1 rad; which is much larger than the 
beaming angle 7-1 < 0.01 (Frail et al. 2001). If we con- 
sider the memory from an axisymmetric jet that is seen 
nearly head-on, not only the memory from the central 
part of the jet vanishes but also the memory added up 
from the edge parts turns out to be canceled because of 
the symmetric distribution of polarization angles. Then 
the GW memory from the acceleration phase of the GRB 
jets can be observed only if the jet is seen off-axis 9.,, > A9 
(Sago et al. 2004). We note that this conclusion depends 
on the geometry of the jet. If the jet is uniform within 
its opening angles, the GW memory from the jet still has 
about half of its maximum values at 9y = A9. 

Similar effects of the finite opening half-angle exist for 
the change in the GW memory from the decelerating 
phase of the GRB jets. For isotropic gamma-ray emission 
in the rest frame of the jet, the GW memory from pho- 
tons is found to be nearly equal with that from the point 
mass that emits the photons ( §2 ). Thus, the memory de- 
creases monotonically owing to the decline in the Lorentz 
factor and the resultant expansion of the anti-beaming 
hole. In this case the GW memory might be observable 
only if 9y > A9 as in the case of the acceleration phase. 

Recent successes in the detection of the polarization 
in the prompt emission of GRBs (Yonetoku et al. 2011) 
suggests the presence of ordered magnetic fields and/or 
an anisotropic photosphere in the emission region. The 
coherence length scale of the magnetic fields is not well 
understood for any GRBs and the origin of such ordered 
magnetic fields is still a debated issue. A variety of differ- 
ent magnetic field configurations and viewing angles have 



B= ^ 

J sin^ 9' cos2 {(f)' -a) + cos^ 9' 

(cos 6*' cos a, cos 0' sina, — sin0'cos(^' ~ a)) . (31) 

We study the case of a = 7r/2, in which the highest degree 
of anisotropy is expected for GW memory. In order to use 
Ah{^,(pB) calculated for the synchrotron model (Equation 
28), we define (I)b such that cosc/)^ = B-t where t is a unit 
vector that is tangent to the great circle connecting the 
points {9',<j)') and {9,,,0) at (6'',0'), 



t = 
where 



p — qcos^ 
sin^ ' 



p = (sin0i,,O,cos6'u), 

q = (sin 9' cos 0' , sin 9' sin (/>' , cos 9' ) , 



(32) 



(33) 



and cos^ = P • q. The angle ^ is equal to 9, the polar angle 
in xyz frame. The GW memory is calculated by summing 
up additive contributions to the amplitude from all points 
in the shell as follows. 



i9'd9'c 



An 



-A/i(e,0s)e 



2irp 



(34) 



where Afi = 27r(l — cos A0) and the integration is per- 
formed over the angles inside the jet. We note that equa- 
tion (12) relates {9,(1}) to (6*', </>'). 

Figure 10 shows the GW memory h = -^/j/i+p + \hx |^ 
calculated for the case of synchrotron emission with or- 
dered magnetic fields with A9 = 0.1 and a = tt/2. We 
show two plots for cases with different viewing angles, 
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Fig. 10. The amplitudes of GW memory h = iy|h.+ p + |/ixP normalized by Ej/d, where Ej is the total energy radiated in 
gamma-rays, are shown for the case of synchrotron emission with ordered magnetic fields with AO = 0.1 and a = w/2. The two panels 
are GW memory for different viewing angles 6^=0 (left) and 0y = 0.001 (right). The time t is arbitrarily scaled. The gamma-ray 
pulses, also plotted with an arbitrary scale, correspond to the collision events, which are modeled by internal shocks with = 100 
relativistic shells of Lorentz factors 7 = 10^ — 10^. 



9v = (left) and 9^ = 0.001 (right). The time is arbitrar- 
ily scaled. The time which passed till the memory reaches 
the plateau corresponds to the duration of the gamma- 
ray emission of the GRB. The amplitude is normalized by 
E^/d, where Ery is the total energy radiated in gamma- 



rays. 



The value of /ijot is shifted so that the final level 



of hjct turns out to be zero. The gamma-ray pulses, also 
plotted with an arbitrary scale, correspond to the collision 
events. We see that the memory from photons emitted at 
each collision is added up, leading to a monotonic rise in 
net memory with time. For 9y=0 the memory from the jet 
is suppressed and unobservable throughout. We ignored 
the broadening of the gamma-ray light curves that should 
appear because of the finite opening angle of the jet. The 
rise of the GW memory at each collision should also be 
smeared. However, even if such smearing is taken into 
account, the waveform, which is made from cumulative 
contributions, would not change its shape so much. 

The behavior of gravitational waveforms shown in 
Figure 10 can be understood as follows. As we changed 
terms in equation (21) as equation (22), we rearrange 
terms in equations (28), (29), and (30). Then, by leav- 
ing only dominant terms, we have 



Ahi9,,cj)B)^^ 
a 



<2{l + cos9)e^"^ > 



2sin^ 



1 — Pm cos 9. 



where E = mri'fr - Im) + msils - 7m)- 
As is found in §2, Ah{9y,4iB) takes a non-zero value of 
order ^ E/d only for ^9^ less than ^ 1 (Figures 7 and 8). 
Thus, we use equation (35) in (34) and change variables 
of integration into 9, (j) to obtain a rough estimate 



(36) 



where E-^ is the energy radiated as photons. This is valid 
both for a single collision and for the sum of radiated 
energy for all the collisions. For = 0.1 and 7 ~ 3 x 
10^ we have h ^ lO^^E^/d. This is a change of the GW 
memory. The initial amplitude of the GW memory is 
carried by the uniform jet and its approximate value is 
given by 



ho . 



'jet 



d A6'2 + 91 ' 



(37) 



where E-^^t is the kinetic energy of the relativistic jet. We 
have confirmed that this approximate expression agrees 
with the exact value (Sago et al. 2004; Hiramatsu ct al. 
2005) to an accuracy of 10 percent for 7 = 10^ — 10"^ and 
a wide range of A0. If we assume Ejct = E^, we obtain 
ha^Ax IQ-^E^/d for = 0.1 and 9^ = 0.001, being in 
a qualitatively good agreement with the results shown in 
Figure 10. 

As is seen in Figures 7 and 8, the change in the GW 
memory is seen only within small angles of 7"^ from the 
line of sight. As long as the small solid angles around the 
line of sight ( a shaded circle in Figure 9 ) is included in 
the jet with an opening half angle A0, which is equiva- 
lent to 9„ < A9 — 1/7, we would observe almost the same 
variation in the GW memory regardless of the position of 
the line of sight. This implies that an ordered magnetic 
field is required only within a small solid angle from which 
gamma-rays we observe are radiated. 
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Fig. 11. The characteristic amplitude hc{f) of the GW 
memory from the decelerating phase of a GRB jet with 
E-y = 10^^ erg. Synchrotron emission from electrons in or- 
dered magnetic fields is assumed. Solid lines represent hc{f) 
for two cases of (d, A6») = (IMpc, 0.1), (lOkpc, 0.1). We as- 
sume a = ■7t/2,9v = for both cases. The amplitude scales 
as A9^^, since we fix Ej and only a small portion of the jet 
of angular size ■y~^ contributes to the change in the GW 
memory. The amplitude does not depend on 6^ as long as 6„ 
is relatively small compared to A6. We have not shown cases 
with dv > AO, since we would then have a significantly small 
variation in the GW memory so that such cases are out of 
interest. 

4. Discussion and Summary 

We studied the GW memory from a radiating and decel- 
erating point mass. We calculated the memory from the 
photons averaged over angular distribution of the emis- 
sion. For isotropic emission (in the rest frame of the 
point mass), the averaged memory has quite similar de- 
pendence on the viewing angle as the memory from the 
point mass. That is, the memory from the photons is 
strongly suppressed for small viewing angles. However 
if the photons are emitted anisotropically, the averaged 
GW memory may have a large amplitude even for small 
viewing angles, enabling us to observe the GW memory of 
GRB jets. We demonstrated an example of the gravita- 
tional waveform expected for the internal shock model in 
the case of synchrotron emission with an ordered magnetic 
field. We find that the GW memory shows a continuous 
rise over the time scale of the gamma-ray emission as well 
as tiny jumps correlated with gamma-ray pulses. Such 
an anisotropic emission of gamma-rays has been inferred 
from the detection of polarizations in the GRB emissions 
(Steele et al. 2009; Yonetoku et al. 2011). Obviously, we 
need more samples of such polarization observations and 
further understanding of the relativistic jet's configuration 
and the mechanism of gamma-ray emission to consider re- 
alistic modeling of the GW memory from GRB jets. 



The GW memory would also provide clues to discrim- 
inating different models of GRBs and related GRB phe- 
nomena. Based on the unified model of GRBs the rela- 
tivistic jet is composed of a number of sub-jets launched 
from a central engine (Yamazaki et al. 2004). In this case 
each sub-jet would be threaded with ordered magnetic 
fields that might be oriented in various directions. Then 
the polarization angles of the GW memory from each sub- 
jet may be different, which leads to a sharp fluctuation of 
the memory with a relatively small amplitude. The unified 
model also predicts that relativistic jets seen with large 
viewing angles are observed as XRFs. For such cases with 
moderate viewing angles the GW memory should vary ow- 
ing to the broadening of the anti-beaming hole. Whether 
the memory has such features in its time variation could 
help us to test GRB models. 

The GW memory from decelerating GRB jets is ex- 
pected to have characteristic frequencies 0.1 — 10 Hz, cor- 
responding to the duration of GRBs, T ^ 0.1 — 10 sec, so 
that it is a suitable target for DECIGO and BBO (Seto, 
Kawamura, & Nakamura 2001; Sago et al. 2004) . Figure 
11 shows the characteristic amplitude hdf) = /|^(/)|, 

h{t)e-'^'''f*dt, from the decelerat- 

-oo 

ing phase of a GRB jet with E-^ ~ 10^^ erg and a dura- 
tion of 1 sec. The noise amplitude hn{f) = [5fSh{f)]^^^, 
where Sh{f) is the spectral density of the strain noise in 
the detector at frequency /, is also shown for Advanced 
LIGO, LISA, and DECIGO/BBO. We use the same for- 
mula for hn{f) or Sh{f) with the one used in Sago et al. 
(2004) and Suwa & Murase (2009). The flat spectrum 
at around / = 1 — 10^ Hz is a direct consequence of the 
time variability of the gamma-ray emission. A GRB with 
= 10^2 _ 10^3 erg and A6' = 0.1 at d = 10 kpc (Galactic 
Center) is likely to provide an amplitude at a level close to 
the sensitivity of Advanced LIGO detector at around 100 
Hz (Flanagan & Hughes 1998). Thus the GW from GRB 
jets may be an interesting target for Einstein Telescope 
(ET) (ET-webpage ) as well as for Advanced LIGO. 

Unfortunately, the local GRB event rate has been 
estimated to be relatively small as ^ 1 Gpc~'^ yr~^ 
(Matsubayashi et al. 2005; Wanderman & Piran 2010). 
This translates into an event rate ^ 10^^ yr^^ ( within 
d=l Mpc ) or 10"^ yr^i ( within d = 10 Mpc ) as a 
rate of GRBs with detectable GW memory. The possi- 
bility of detecting the GW memory from on-axis viewing 
angles, that is, almost simultaneously with gamma-rays 
from GRBs, makes it easier to extract the memory com- 
ponent embedded from the detector's signals. Thus the 
GW memory from GRB jets will be one of the interesting 
targets for the next generation detectors such as KAGRA 
(Somiya et al. 2012), DECIGO/BBO, and ET. 
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